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ABSTRACT 

The growth of galaxies through adiabatic accretion of dark matter is one of the main drivers of 
galaxy evolution. By isolating it from other processes like mergers, we analyse how it affects 
the evolution of star clusters. Our study comprises a fast and approximate exploration of the 
orbital and intrinsic cluster parameter space, and more detailed monitoring of their evolution, 
through iV-body simulations for a handful of cases. We find that the properties of present- 
day star clusters and their tidal tails differ very little, whether the clusters are embedded in a 
growing galactic halo for 12 Gyr, or in a static one. 
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1 INTRODUCTION 


Since the time of their formation, galaxies have undergone a 
variety of transformations, from major mergers to slow accre¬ 
tion of dark matter and intergalactic gas l Press^^^cheghtejQ^TJ; 
IWhite & Reeslll978l : Ike Fevre et al.ll2000l: IPekel et al.ll2009l) . As 
old objects evolving in these galaxies, globular clusters and dwarf 
galaxies can pote ntially probe such changes and record some 
of the ir signatures dKrauss & Chabov^l2003l : lKravtsov & GnedinI 
l2005h . It is thus important to understand how the evolution of galax¬ 
ies has impacted present-day stellar populations. Ideally, theoretical 
studies on this topic should account for the coupling between the 
internal dynamics of small stellar systems (clusters and dwarfs), 
and the effect of their galactic environment in cosmological con¬ 
text. However, because of the wide range of scales and of physical 
processes involved, considering hoth aspects simultaneously can be 
challenging. 

One possible approach consists in focussing on low-density 
objects like dwarf galaxies, where star-star interactions can be ne¬ 
glected. By adopting a collisionless treatment of gravitation in nu¬ 
merical simulations, one can efficiently explore a wide parameter 
space. Such methodology is la rgely used in the studies of satellite 

galaxies and tidal streams (e.g. Penarrubia et al .11 20051: IBov 


iBonaca et a Dl2014l:lKuDPer et alj2015l)~ 

But for denser systems like star clusters, two-body relaxation 
plays a paramount role in setting the flux of st ars escaping the clus¬ 
ter, a nd thus in the evolution of the cluster jpukushige & Heggiel 
bOOOh . Therefore, gravitation must be treated in a collisional fash¬ 
ion, which involves an important numerical cost, a nd forbids to d o 
so over galactic and cosmological scales (see also lDehnenll2014h . 
Therefore, the large scale influence is not (yet) considered self- 
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consistently in star cluster A^-body simulations. Up to now, several 
paths have been followed to implement the tidal effects from the 
galaxy on a cluster: 


(i) the cluster remains within a single galaxy for its en- 


tire lifetime (e.g 

[Baumgardt & Makinoll2003l: 

Kunner et al.l 

201 c; 

Hurley & Sharal 

2 OI 2 I; Madrid et all 1201 2l; Vesnerini et all 

2014; 

Webb et alj2014 

). The galactic potential is usually static, and often 


axisymmetric. 


(ii) the cluster is accreted from a satellite galaxy onto a mas¬ 
sive galaxy. The accretion event is often mod elled by replacing one 
galactic potential wi th another, both static dMiholics et aDl2014l : 
ISianchini et ahlFoiSl) . 

(iii) the galactic potential is evolved using a separate galaxy sim¬ 
ulation to allow for time-evolving, n on-analytical potentials and 
include effects like galaxy mergers l lRenaud & Gielej l2013h or 
even complex tidal histories in cosmological context jRieder et al.l 
|2013l) . 

Although the last approach provides a relatively high level of 
realism, it does not allow us to distinguish the relative role of the 
numerous physical processes involved in the evolution of the clus¬ 
ters. In this paper, we focus on a specific aspect of the galactic tides, 
namely the secular, adiabatic cosmological growth of galaxies, and 
neglect other effects, at the expense of realism. The full story of the 
co-evolution of star clusters and their hosts would be told by con¬ 
sidering a combination of this particular effect, and the other mech¬ 
anisms driving galaxy evolution, like minor and major mergers, and 
the accretion of gas. In that respect, cosmological simulations and 
merger trees would serve as a base to establish the relative weight 
and frequency of these events, and to allow us to build the true tidal 
history of clusters, in future studies. 

Here, we follow the evolution of clusters in a time-dependent 
galactic potential, and compare to that in static potentials, to quan- 
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tify the role of secular galaxy evolution on the present-day proper¬ 
ties of star clusters. Our method is two-fold: (/) an exploration of 
the parameter space performed with very fast codes, hut at the price 
of some simplifications, and (h) several much slower hut more ac¬ 
curate A'^-body simulations using the relevant sets of parameters 
identified in the first step. 


2 METHODOLOGY 
2.1 Time-evolving potential 


For simplicity and to limit degeneracy in the results, we only con¬ 
sider the halo component of the galaxy and opt for a self-similar 
growth. The galaxy does not experience any merger event but 
slowly and smoothly grows (in mass and radius) with time. Fur¬ 
thermore, the halo remains spherically symmetric and we neglect 
dynamical friction. We choose to make such simplifications in or¬ 
der to focus on the role of a single physical process (namely the 
adiabatic cosmological growth). 

We use the analytica l description o f a growing 
iNavano. Frenk & Whitel ( Il997li halo proposed by ISuist & Helmil 
(l2014h , who fitted t he evolution of ha lo parameters using the 
Aquarius simulation dSpringel et alj2008li . Namely, the mass-scale 
and scale-length of the halo evolve with redshift a as 

f Ms[z) = Ms,o exp (-0.2z) 

\ Rs{z) = i?s,o exp (—O.la). 


We have adopted the values of 0.1 for th e mass growth parame - 
ter (flg) and of 2.0 for the 7 parameter in iBuist & Helmil ^2014) . 
their equations 22 and 23. This provides an evolutio n comparable to 
that of the Milky Way-like halo labelled Aq-E-4 in iBuist & Helmil 
( |2014 their Fig. 4). The effect of changing these parameters is dis¬ 
cussed in Section[5] The galactic potential, as a function of radius 
r and redshift « is then 


4'a{r,z) 


GMsjz) 

r 


In 



( 2 ) 


(where G is the universal constant of gravitation), and corresponds 
to the density profile: 


Poir, z) 


Ms{z) 

dvrr [r -f Rs{z)f 


(3) 


To evaluate the paramaters of the potential during the simula¬ 
tion of the cluster, we convert the time t since the beginning of the 
simulation (i.e. the age of cluster) into redshift using 


1 + z 


f i-n, 

\ Gjti 






(which follows from equation 13.20 of |Peeblesll993ll . to represents 
the age of the Universe when the simulation is started. We adopt 
the values of the co smological constants firn ~ 0-3 1 and Ho — 
68 km Mpc”*^ jPlanck Collaboration et al.l2014h such that the 
Hubble time equals 13.7 Gyr. 

We start our study at redshift 5, which leads to to ~ 
1169 Myr, and we consider the evolution of clusters over 
12.56 Gyr (corresponding to the time difference between z = 5 
and z — 0. We choose the present day (z — 0) values of the 
mass-scale and scale-length to be Ms,o = 1.5 x 10^*^ M© and 
Rs,o = 16 kpc. Fig.[T]shows the evolution of the mass-scale (Ms) 
and scale-length (Rs), both normalised to their values at z = 0. 


5.0 3.0 2.0 1.5 1.0 ^ 0.5 0.3 0.1 0.0 



Figure 1. Evolution of the scale parameters of the galaxy, normalised to 
their value at z = 0 (i.e. for t Ki 12.6 Gyr). t = 0 corresponds to the start 
of the simulations (i.e. z = 5). 


2.2 Galaxies 

Our study of the effects of time-varying potentials on the evolution 
of star cluster comprises three types of tidal histories: 

• S5 (static, z = 5): the cluster evolves in a static version of the 
potential of Equation l|3, with properties corresponding to z = 5. 
Equation 0 provides the corresponding galactic mass scale and 
scale radius: « 5.5 x 10'^° Mq and Ri 9.7 kpc. 

• TD (time dependent): the initial setup is strictly the same as 
for S5, but the potential evolves with redshift from z = 5 to z = 0, 
as described in the previous Section. 

• SO (static, z = 0): the cluster reaches the exact same final 
orbital position and velocity than in TD, but has evolved in a static 
potential with properties for z = 0. In practice, we compute the 
final position and velocity of the cluster from TD, when z = 0 
has been reached. We then “freeze” the potential, and perform a 
backward integration of the cluster orbit. This gives us the initial 
position and velocity of the cluster for the SO run. 

Therefore, by setting one initial position and velocity for a 
cluster in the S5 galaxy, we uniquely define three orbits and the 
associated three tidal histories: S5, TD and SO. 

For simplicity, we always initially set the cluster at apocen- 
ter, with a purely tangential velocity, in our S5 and TD cases. 
Orbit integration is pe rformed using the NBODY6tt method 
( iRenaud & Gielesll20151) . either to integrate the motion of the clus¬ 
ter (Section 0, or in the full W-body context where the cluster is 
described star-by-star (Section|4l(. 


2.3 Star cluster fiducial initial conditions 

We considered several intrinsic and orbital initial conditions for 
the clusters. Unless otherwise mentioned, our clusters are modeled 
with 32768 stars distributed on a Plummer sphere with a virial ra¬ 
dius of 3 pc (i.e. a half-mass ra dius of « 2.3 pc). The masses of 
the stars follow a lKrounal l l200lh initial mass function, from 0.1 to 
1 M 0 , leading to a total initial mass of Ri 1.03 x 10“* M©. We do 
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not account for stellar evolution. We explore other cases by varying 
these parameters in the next Sections. 


3 PARAMETER SPACE EXPLORATION 

We aim to determine if, under which conditions and to what ex¬ 
tend, the evolution of a cluster in a cosmologically growing poten¬ 
tial (TD) differs from that of the same cluster in a static z = 0 ver¬ 
sion of this potential (SO). This difference relies on (/) the strength 
of the tidal field and (it) how sensitive the cluster is to tides. 


3.1 Quantifying the tidal field 


We first seek an estimate of the galactic disruptive effect on star 
clusters. Such effect comprises contributions of gravitational (iner¬ 
tial) and orbital origin (non-inertial). For non-circular orbits, non- 
inertial effects, like the centrifugal force, do not yield an analytical 
expression. However, we can estimate their contributions by con¬ 
sidering that the orbit is instantaneously circular, i.e. by neglecting 
the Coriolis and Euler effects, and by computing the centrifugal ef¬ 
fect using the inst antaneous orb i tal ang ular frequency. In this ide¬ 
alised framework. iRenaud et al.l l l201lh provides the expression of 
the effective tidal tensor, which encompasses both the gravitational 
and centrifugal effects. Using the expression of the galactic poten¬ 
tial (Equation^, the main eigenvalue of the effectiv^ tidal tensor 
reads 


Ae(r,z) 


GMsjz) 

rp3 




4r^ -I- SrRs{z) 
[r + Rs{z)f 



This quantity represents the galactic effect along the galaxy-cluster 
axis, and can be used to estimate the tidal radiufl (see lRenaud et all 
IroTU for details). 

We show an example of evolution of the effective eigenvalue 
in Fig. 12 together with the galactocentric radius and the galactic 
mass enclosed in this radius, for the three orbits (TD, S5 and SO). 
As constructed, the TD case evolves from the S5 to the SO setup. 
On average, despite of a smaller orbital radius, the mass enclosed 
within the galactocentric radius in SO is larger than that in S5, be¬ 
cause the SO galaxy yields a higher total mass. We estimate the 
mean tidal strength over an orbit by computing the time-average 
(Ae) of the main effective eigenvalue over an orbital period (dashed 
lines on Fig. |2j- This quantity is constant in static potentials. The 
larger enclosed mass and smaller galactocentric distance result in a 
stronger tidal effect along the SO orbit than in the S5 case. 

The secular cosmological galactic growth can only affect the 
evolution of star clusters if the tidal effects in the S5 and SO case 



20 40 60 80 100 

i's 5 (t = 0) [kpc] 


Figure 3. Relative strength of the tidal field between the SO and S5 ca.ses, 
as a function of initial galactocentric radius {rss{t = 0)) and orbital ec¬ 
centricity in the S5 case (see text for details). The difference between tidal 
strengths is maximal for circular orbits in the outer region of the halo. 


are significantly different, or in other words, if the ratio 

(Ae.So) 

(Ac.ss) 


( 6 ) 


is large with respect to unity. This ratio depends on the initial posi¬ 
tion and velocity of the cluster in a non-trivial way. By varying the 
initial galactocentric radius and orbital eccentricity of the cluster 
for the S5 orbit and integrating it numerically, we obtain the map 
of the ratio of tidal strengths shown in Fig. |2 (Recall from Sec- 
tion l2.2l that setting the initial position and velocity for S5 uniquely 
defines the initial position and velocity for SO.) 

The largest differences between the tidal fields of the redshift 
5 and 0 galaxies are found in the outer regions of the galactic 
halo, and for circular orbits. Secular galactic growth induces the 
largest differences in the outer regions of galactic halos and thus, 
the largest differences in tides are found along orbits that remain in 
such regions for the largest fraction of their period. 

However, at large galactocentric distance, tidal forces are 
weak and only have a mild impact on the evolution of clusters. 
The differences found between the S5 and SO cases might thus not 
translate into differences in the properties of the clusters. 

To summarise, we have identified the orbits favouring the 
largest differences in tides between high and low redshift, but the 
resistance of clusters to tidal harassment must also be considered 
before concluding on the effect of secular galactic growth on star 
clusters. 


^ The main eigenvalue of the purely gravitational tidal tensor (i.e. neglect¬ 
ing the centrifugal effect) is: 


X{r, z) 


GMs(z) 




+ 2rRs(z)\ 

[r + Rs{z)f ]■ 


^ In the textbook context of circular orbits around point-mass es, us i ng the 
ef fective eigenvalue lea d s to th e definition of the tidal radius of lKind H962h 
or iBinnev & Tremaine! i2008l) , while t he pu rely gravitational eigenvalue 
corresponds to the definition of lSnitzeil Il987h . 


3.2 Cluster sensitivity to tidal harassment 

We have seen in the previous Section that the average strength of 
the tidal field experienced by a cluster could increase by a factor 
of a few because of the secular cosmological growth of the galaxy. 
We focus here on star cluster responses to such differences. 

One of the most accurate manner to study star cluster evo¬ 
lution is through A^-body simulations. These are however very 
numerically costly and such approach forbids an exploration of 
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time [Gyr] 


Figure 2. Top: galactocentric radius for the three orbits (S5 in red, TD in black and SO in blue). For the S5 and TD runs, the cluster is initially set 20 kpc 
from the galactic center, with a purely tangential velocity of 50 km s“^. The initial conditions for the SO runs are constructed as explained in Section 
The dotted lines indicate the coiTesponding scale radius {Rs) of the galactic potential. Middle: galactic mass enclosed within the instantaneous orbital radius. 
Bottom: main eigenvalue of the effective tidal tensor, under the approximation that the orbit is instantaneously circular (see text, EquationO. The dashed lines 
indicate the time-average of the eigenvalue along the S5 and SO orbits. 


a wide parameter space. At first order however, cluster evo¬ 
lution re lies on an handful of coupled differential equations 
Ambartsum ani 19381: Chandrasekhaj 19421: Kindll95^: 


:lHeggielll97i iHut et alJ[I^:lLee & Ostrikei||l987|: 
20nh. The code EMAC SS jAlexander & GielMboi/ ' 
I 2 OI 4 I: [Alexander et alJ|20lT solves these equations 
numerically and provides an easy and very fast way to evaluate 
the properties of the cluster along its evolution. 



We first setup our fiducial cluster (Section [2.3b in EMACSS, 
using the tidal field strengths computed in the previous Section. 
To model the tidal effect, we determine the mass of the point- 
mass galaxy that would lead to the same average effective eigen¬ 
value (Ae) than our NFW halos, at a given orbital radius. Do¬ 
ing this allows us to set the same tidal radii in EMACSS than the 
time-average ones measured in the previous Section. However, we 
neglect the differences betw een the shape of the potentials (see 
iTanikawa & FukushigdboiT for a discussion on that matter). This 
assumption will be validated later. 


Fig. a shows the ratio between the final cluster mass when 
embedded in the SO and S5 tides. Despite the differences in tidal 
strengths found in Fig. [3] the differences in the final mass of the 
cluster remain below 10 percent in all cases. In other words, the 
final mass of our fiducial cluster is almost the same in S5 and SO 
for all initial orbital radii and eccentricities considered. 


To extend this conclusion, we apply the same method to clus¬ 
ters with different sensitivity to tides by varying the initial half¬ 
mass radius and thus, indirectly, the Roche-filling factor (the ratio 
between the cluster half-mass radius and the tidal radius). We plot 
in Fig.|5]the ratios of the final masses (as in Fig. lU and final half¬ 
mass radii. We arbitrarily choose the tidal field strength of orbits 
with an eccentric of 0.3, but reach the same conclusion of other 



Figure 4. Ratio of the final cluster mass (computed with EMACSS, see text) 
when evolved in the SO and S5 potentials, as a function of initial galac¬ 
tocentric radius and orbital eccenh'icity in the S5 case. The cluster initial 
parameters are those of the fiducial case. Note that EMACSS is not designed 
to treat orbits with a high Roche filling factor, and thus the cases within the 
central 5 kpc are not considered here. 


values, as already suggested by the independence to eccentricity 
showed in Fig.|4] In all cases explored, the properties of the cluster 
(mass, size) along the S5 orbit lie within less than a few percent 
from those along the SO orbit. 

Because the TD case represents a slow and smooth transition 
between S5 and SO, the differences between the latter can be seen 
as an upper limit of the expected differences between the time- 
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Figure 6. Evolution of the mass (top) and the Lagrange radii (10, 50 and 
90 percent of the mass) of our fiducial cluster along the orbits showed 
in Fig. [ 2 ] computed with NBODY6tt. Only stars with a negative energy 
(kinetic + potential from other cluster members) ar considered here, as in 
iRenaud & Gielej2015h . 


Figure 5. Top: ratio of the final cluster mass when evolved in the SO and S5 
potentials, as a function of initial orbital radius (with an eccentricity of 0.3) 
and initial half-mass radius. The black symbols correspond to the 7V-body 
simulations described in Section [ 4 ] Bottom: same but for the final half¬ 
mass radius. The horizontal ridge visible in this plot is an effect of clusters 
reaching the core-collapse phase shortly before the end of the simulation 
(5=iJ 12.6 Gyr). Clusters with a long relaxation time (large half-mass radius, 
above the ridge) are still in the unbalanced pre-core-collapse phase, with 
a de creasing half-mass radius, when the EMACSS simulations are stopped 
(see lGieles et alj20l4 for details). 

dependent (TD) and static (SO) cases. Since this upper limit is very 
small, we expect clusters with a time-dependent tidal history to 
share very similar properties as those in static potentials. In other 
words, cluster evolution should be fairly independent of the evolu¬ 
tion of its host galaxy (for the secular growth considered here). 


4 A^-BODY SIMULATIONS 

When performing the exploration of the parameter space described 
in the previous Section, we made the assumption that the time- 
average tidal strength can be used to infer the galactic influence on 
star clusters. While this is perfectly valid for circular orbits, it is no t 
for eccentric cases jBaumgardt & Makinol2003l : IWebb et alj2014h . 
An analytical description of tidal effects for non-circular orbit is yet 
to be established. Furthermore, we have neglected the second order 
effect of the detailed shape of the tidal field when assuming the 
tidal radius is a good representation of the tidal strength. 

To test these assumptions and validate our conclusions, we 
select a few cases from our parameter space study, and perform 


full A^-body simulations. We use the method and im plementa- 
tion NBO DY6tt presented iniR enaud & Gielej (f2015h . based on 
NBODY 6 l lAarsethl2003l : lNitadori & AarsethI20121) . The method re¬ 
lies on a description of the galactic potential (through a numerical 
routine) as a function of position and time. Using this definition, 
the code integrates the motion of the cluster around the galaxy, and 
adds tidal acceleration to the internal acceleration for all stars. We 
have defined the galactic potentials and their time-dependence ex¬ 
actly as describe in Section ITT] and performed several sets of S5, 
TD and SO runs. 


4.1 Cluster structure 

Fig.|6] shows the evolution of the mass and some Lagrange radii 
of our fiducial cluster (Section 12.3b set at an initial position of 
20 kpc with an orbital eccentricity of 0.3 (as in Fig.^. As expected 
from the parameter space exploration presented above (black dot 
in Fig. O, the relative difference in final mass (respectively half¬ 
mass radius) between the S5 and SO runs is very small (« 6.1 per¬ 
cent, respectively ~ 4 percent). The excellent agreement between 
the NBODY6tt results and EMACSS validates our assumptions, at 
least in this region of the parameter space. We have also tested this 
agreement by considering other eccentricities (0, 0.5 and 0.7, not 
shown here) and reached the same conclusions. As foreseen in the 
previous Section, the difference between the SO and TD runs is 
even smaller (« 2.6 percent for the mass and ~ 3 percent for the 
half-mass radius) than that between the S5 and SO cases. 

We also run other NB0DY6tt models, corresponding to the 
black triangles on Fig. and reach the same conclusions, both on 
the validity of our method, and on the physical results obtained. 
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Figure 7. Left: Position of the stars at redshift 0 for the TD (black) and 
SO (blue) runs. The plus sign indicates the position of the center of the 
galactic halo. Right: Isodensity contours of the clusters and the tidal tails. 
The dashed line represents the orbit of the cluster. The two models overlap 
almost perfectly. 



galactic azimuth [°] 


Figure 8. Dishibution of mean galactocentric distance of stars (top) and of 
the mass (bottom), as a function of azimuth in the galactic reference frame, 
for clusters having evolved s; 12.6 Gyr in the TD and SO potentials. (The 
cluster are centered on the azimuth 0°.) The dashed line represents the orbit 
of the cluster. 


4.2 Tidal tails 

Tidal ejection of stars from clusters leads to the formation of tidal 
tails. Although we have found that the mass loss rate of clusters 
appear to be independent of tidal history (in the context of secular 
growth), tidal tails could respond differently to a time-dependent 
and a static galactic potential. 

Fig.|7]shows the position of all the stars (cluster and tails), and 


Fig.[8]displays the distribution of mean galactocentric distance and 
the mass, in galatic azimuth bins (of width 1°), for the clusters hav¬ 
ing experienced the TD and SO tidal histories. Differences between 
the models appear at rather large distance from the cluster center. 
The differences in mass are of the order of a few solar masses, i.e. 
concern only an handful of stars. 

Despite different tidal histories over the 12.6 Gyr of evolution 
we considered, the final position, density, length and distribution of 
substructures in the tidal tails between the TD and SO models are 
almost undistinguishable. In other words, clusters and their tidal 
debris do not retain signatures of the tidal history they experienced. 


5 DISCUSSION AND CONCLUSIONS 


We study the evolution of star clusters embedded in cosmologically 
growing galaxies, and focus on the secular, adiabatic growth of dark 
matter halos, neglecting impulsive and transient events like galaxy 
mergers. Our main findings are: 


• Although the typical tidal fields associated with high and low 
redshift halos can vary by a factor of several, these variations do 
not translate into major differences in the properties of star clusters 
embedded in those fields. 

• Because of these similarities, the details on the secular evolu¬ 
tion of the halo (growth rate, growth epoch) have a negligible effect 
on the present-day properties of the clusters. 

• Present-day star clusters that co-evolved with their host galaxy 
yield the same properties as if they had evolved in a static halo. 


Clusters orbiting in the innermost regions of galaxies are af¬ 
fected by the tidal effect from the baryonic components that we 
have neglected. Among other effects, the details on the formation 
of thin and thick discs could modify the role of disc shocking in the 
evolution of clusters ( ISpitzeiiri958h and alter our conclusions. We 
still understand too little on galaxy and structure formation to reach 
conclusion on this matter. 

Furthermore, we only considered spherically symmetric ha¬ 
los, i.e. neglecting (time-dependent) anisotropy and substructures. 
Although we expect t hese aspects to rarely be of significant impor¬ 
tance for star clusters. iBonaca et al.l ( l2014l) showed that could alter 
the morphology of tidal debris. 

Our conclusions depend on the growth history of the galaxy, 
mainly when and how fast the bulk of the growth takes place. Cos¬ 
mological simulations indicate that most o f the adiabatic growt h 
happens at high redshift (recall Fig. [T] see iBuist & FIelmill2Q14l) . 
which has a very mild effect on star clusters, as we have shown. A 
galaxy can also experience a significant growth in the form of a vio¬ 
lent event. If so, the growth enters the impulsive regime and can be 
seen a galaxy merger. The response of th e clusters to such event i s 
non-negligible but complex, as showed in iRenaud & Giel^ ( l2013[) . 
The real evolution of a galaxy and of its clusters lies in between 
these two extremes. Our study demonstrates that the main impacts 
of galactic tides on star clusters in the cosmological context are the 
impulsive events, and not the adiabatic growth. 

In the context of the Milky-Way, the abse nce of evidence for 
recent major mergers (in the last ~ 6 — 9 Gvr.lP eason et ^ J2013l 
and/or since the formation of the disc iRuchti et al.ll2014 ) makes 
our results directly applicable to clusters formed in-situ. Studies of 
such clusters can thus safely make the assumption that the galac- 


(© 0000 RAS, MNRAS 000, 000-000 




























































Secular galactic growth on star clusters 1 


tic potential has been static over several Gyr, without biasing the 
conclusions on the mass- and size-functions of the clusters. Stud¬ 
ies that consider globular cluster population evolution in a static 
Milky Way potential showed that it is not possible, within a Hubble 
time, to evolve an initial cluster mass function with a power-law 
shape and an index of -2 (as found for young massive cluster today, 
[Portegies Zwart et alj201oll . into a mass function that is peaked at a 
universal mass of ~ 2 x 10^ Mf?) with secular dynamical evol ution 

( lBaumgardlll998l : IVesDerinj200lMGieles & Baumgardtl2008l) . Our 

results support this conclusion. Note however that this is not valid 
when considering the non-negligible fraction of clusters of external 
origin, which h ave been accreted from dwarf satellite galaxies onto 
the Milky Way jMarm-Franch et alj2009l : lLeaman et alj|2013h . 

The diversity of scenarios of galaxy evolution seen in cos¬ 
mological simulations implies that present-day star clusters have 
experienced a wide variety of tidal histories. By decomposing 
these scena rios into individual pro cesses and events like ma- 
jor merge rs dRen aud & Gielesll2013h . accretion of dwarf satellites 
dBianchini et alj|2015l) . secular growth (this work) and others, and 
by understanding their relative roles on the evolution of clusters, 
we will soon able to seek and identify specific imprints of galaxy 
evolution on the observed properties of star clusters. 
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